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Assuming an adiabatic evolution of a gamma-ray burst (GRB) fireball interacting witli 
an external medium, we calculate the hydrodynamics of the fireball with energy injection 
from a strongly magnetic millisecond pulsar through magnetic dipole radiation, and obtain 
the light curve of the optical afterglow from the fireball by synchrotron radiation. Results 
are given both for a homogeneous external medium and for a wind ejected by GRB 
progenitor. Our calculations are also available in both ultra-relativistic and non-relativistic 
phases. Furthermore, the observed R-band light curve of GRB000301C can be well fitted in 
our model, which might provide a probe of the properties of GRB progenitors. 

PACS: 98.70.Rz, 97.60. Jd, 95.30.Lz 

Since the measurement of the redshift of GRB970508, gamma-ray bursts(GRBs) as the 
cosmological origin are confirmed. '^1 Thus, one of the most important issues is the nature 
of the central engine. The popular theoretical explanation of their radiation properties 
has been commonly thought to be the fireball model t^^, in which GRBs result from the 
dissipation of the kinetic energy of the relativistically expanding fireballs. Two main ideas 
have been proposed to realize the dissipation: the internal and external shock model. In 
present letter, we will consider the external shock model which reproduces very well the 
delayed emission at lower energy. 

A number of afterglows from GRBs have been observed at X-ray, optical and radio 
wavelengths, whose major features can be explained by simple standard models. '^''^1 We 
consider a newborn millisecond pulsar with a strong magnetic field at the center of the 
fireball, which may be formed by several models: accretion-induced collapse of magnetized 
white dwarfs, merge of two neutron stars, and accretion-induced phase transitions of 
neutron stars to strange stars. '^'^1 It has been proposed that very strongly magnetized 
pulsars may be the central engine of GRBs.^^^ It is natural to expect that if a GRB results 
from the birth of a strongly magnetic millisecond pulsar, then after the main GRB, the 
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pulsar continuously supplies energy to the fireball through magnetic dipole radiation. And 
during the process, we haven't consider gravitational radiation due to the very high surface 
temperature of the newborn pulsar. The power of the pulsar is radiated away mainly 
through electromagnetic waves with frequency of a; = 27r/P, where P is the period of 
the pulsar. Once the electromagnetic waves propagate in the shocked interstellar medium 
(ISM), they will be absorbed for ujp > a;,^^^ where cUp is the plasma frequency of the shocked 
ISM. 

At the center of the fireball, the pulsar loses its rotational energy through magnetic 
dipole radiation, [''1 whose power is given by 



L = 4 X lO''Bl,,Pr^,R%s,, + erg s-\ 
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where B±^i2 — BsSinO/lO^Hn units of G, Eg is the surface dipole field strength, 9 the angle 
between the magnetic dipole moment and rotational axis, Pi^ms the initial period in units of 
1ms, and Rns,6 is the neutron star(NS) radius in units of lO^cm, tj, is the burster time and 
T is the initial spin-down timescale defined by T = P/2P. 

In the following we will consider the adiabatic expansion of a fireball. In the fixed 
frame, the total kinetic energy of the fireball is'^^l 

Ek={-i- l)(Mo + M,^y + 7t/, (2) 

where 7 is the bulk Lorentz factor, Mq is the rest mass ejected from the GRB central 
engine, Mgw is the rest mass of the swept-up medium, U is the internal energy. As Huang et 
g^|ji2,i3] suggested f/ = (7 — \)Msw(?-, available in both ultra-relativistic and non-relativistic 
phases. For an adiabatic fireball, the change of the kinetic energy equals to the energy 
which the fireball has obtained from the pulsar through magnetic dipole radiation: 

J 771 

^ = (1-W), (3) 

where (5 — (1 — 1/7^)^, t — tj, — R/c is the observed time, and R the blast wave radius. 
Prom Eqs (2) and (3), we can find 

dr^ ^ {l-f3)L{t) _ 7' - 1 dMs^ 

dtb (Mo + 27M,^)c2 Mo + 27M,^ dh ' ^ ^ 



Generally, it is assumed that the medium surrounding the GRB source is homogeneous, 
n{R) — n*. However, if a collapsing massive start-*^^! is the origin of a relativistic fireball, the 
circum-burst medium is the wind ejected by the star prior to its collapse, whose density 
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decreases outwards. Recently, the discovery of the connection between supernovae and 
GRBs is a strong support to the assumption that some GRBs come from coUapsars.'^^l 
The dynamical evolution and afterglow in the wind model are expected to be very different 
from those in the homogeneous case. Significant work in this direction has been done by 
some researchers, f^'^^] The baryon number density of the wind medium is expected to be 
n{R) = AR-"^, A is scaled to A = 3.0 x lO^^A^cm-^I^^l and A^ = M/v where M is the 
mass loss rate of massive star in unit of MQyr"^ that has ejected the wind at constant 
speed V in units of lO^kms"^. In present work, we give more detailed calculations in both 
the homogeneous external medium and the wind, and investigate the difference between the 
dynamical evolution and light curves of afterglows arising for two types of external medium 
to find ways for distinguishing the two models. 

Based on the model described above, we have evaluated the propagation of the blast 
wave numerically by taking Eq — lO^-'^ergs, = lcm~^, A^ — 1, Mq = 5 x 10~^Mq, 
Pi — 1ms, Ens — lO^cm, and 7 = 2.0 x lO'^^g cm^. To describe completely the properties 
of the case with energy injection, we also calculate afterglows without energy injection 
as a comparison. In Figs. 1 and 2, the solid lines correspond to the case without energy 
injection, the dotted lines to the surface magnetic field strength of the pulsar, Bg = lO^^G, 
the dashed lines to Bg = lO^^G and the dash-dotted lines to Bg = lO^'^G. 

Figure 1 shows the evolution of the bulk Lorentz factor in the homogeneous ISM and 
wind models. From this figure, we can see that the bulk Lorentz factor decreases faster 
firstly but more slowly later in the wind model than in the homogeneous case. Comparing 
the cases of with and without energy injection, we find that the bulk Lorentz factor in 
the former is slightly larger than that of the latter. With different surface magnetic field 
strengths of the pulsars, the timescales at which the effect of energy injection on the 
evolution of the fireball becomes significant vary accordingly. 

To compare with observations, we need calculate synchrotron radiation from the 
shocked medium, in which we take the same electron distribution as in Dai et aL.'^^l In 
Fig. 2, we show the computed R-band optical light curves. In these figures, some of the 
parameters are fixed (p = 2.5, = 0.1, = 0.01, — 3Gpc). Figure 2a corresponds to 
the homogeneous case, and Figure 2b to the wind case. We clearly see that the afterglow 
light curves in the wind model are steeper than those in the homogenous medium, verifying 
our expectation. Furthermore, we find that the magnitude of the optical afterglow from 
the fireball first decreases rapidly with time, subsequently fiattens, and finally declines 
again. From the existence of the fiattening in the light curves, we may successfully explain 
the flattening and steepening optical features observed in some GRB afterglows which 
contradict the results of the standard model. 



-4- 



According to the numerical results above, we can use our model to study some 
GRBs whose afterglows appear anomalous, such as GRB 000301C. We employ a standard 
cosmology with Ho = 65km s-^Mpc-^, = 1-0 and A = 0. GRB 000301C is one of the 
latest afterglows exhibiting a sharp break in the optical light curve. The redshift was 
measured using the Hubble space telescope to be 1.95 ± 0.1 by Smette et al.'^^^ and was 
later refined by Castro ct al.'^^l using the Keck 10-m telescope to a value of 2.0335 ± 0.0003, 
corresponding to Dl ^ 13Gpc. Figure 3 shows our numerical fit to the R-band light curve of 
GRB 000301C. We note that the radio afterglow of GRB 000301C shows a break around 10 
days after the burst, which is similar to the radio afterglow of GRB 980519. This suggests 
that the circum-burst environment may be a dense medium. l^^l Dai et al.^^'^^l have studied 
the dense medium model of GRB afterglows in detail. Hence we use our homogeneous 
model with the following fixed values: = 1.0 x lO^^ergs, Mq = 5.0 x 10~^Mq, Pi = 
0.5ms, = 5 X lO^cm-^, p = 3.0, = 0.02, = 0.0005, = 3.0 x IO^^q ^nd Dl = 13Gpc. 
The R-band afterglow can be described by a power-law with the index a ~ 1.2 when 
1.3 < t < 7.3d and a ~ 2.7 when 7.3 < t < 31. 5d, and the break at t = 7.3d which is well 
fitted in our model. 

In this letter, the significant results of our calculation show that afterglow light curves 
may have two breaks. Then our model should be well fitted to the light curves of GRB 
afterglows exhibiting one break or two breaks, which depart from the explanation of the 
standard model. Recently, Dai & Lu'^^lhave used this model to analyze the unusual optical 
afterglow of GRB 000301C. The first break occurs at the time when the energy injection 
becomes significant, and this timescale is estimated by r 5 x 10'^ Eo^5iBj^^2Pi,ms^6^^- 
When we take different magnetic field strength {Bs = 10^^, 10^^ and lO^^G), r corresponds 
to the different timescales (r = 5 x 10^, 5 x 10^ and 5 x lO^s) which are consistent with our 
calculated results. 

We have numerically compared two types of external medium: homogeneous and 
wind cases. Note that the light curves in the wind model are steeper than those in 
the homogeneous case, so one may use the wind model to explain some afterglows that 
have decayed rapidly at later times which may provide a way to distinguish the two 
cases, giving some hints in understanding GRB progenitor environments. However, the 
sideways-expansion effect of a jet can lead to similar steep light curves, so one 
urgently know how to further distinguish the wind and jet effects. Now it is widely argued 
that the break may be due to a jet-like outfiow at 7 ~ 9~^, where 9 is the half opening 
angle of the jet. 

During our numerical calculations, if we enlarge the value of initial fireball energy, we 
find that the light curves have no apparent fiattening even we change the surface magnetic 
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field strength of a pulsar. We estimate the total energy injection from the pulsar through 
magnetic dipole radiation, Etot = 2.0 x lO^^P'j^^I^^ in units of erg , which is independent of 
the magnetic field strength. When we take Pj = 1ms and / = lO^^g cm^, Efot is only about 
lO^^erg, which is negligible compared to the initial fireball energy up to 10^^ or larger. Since 
the moment of inertia is in a narrow range for different equations of state (EOS) at high 
density, we have to decrease the period of the pulsar if we need to explain the afterglows 
of some GRBs with high initial energy observed. Note that the submillisecond period for 
a neutron star is difficult to achieve using reasonable EOSs for normal neutron matter. I^'^l 
However, strange stars can reach much shorter periods than neutron stars due to the 
existence of high viscosity in strange stars. [^^1 If the pulsar with submillisecond period would 
exist, it may be a strange star rather than a neutron star. We fit the R-band light curve 
of GRB 000301C with Pi = 0.5ms, implying that the pulsar, if it was the center engine of 
the fireball, may be a strange star in our model. Therefore, some unusual afterglows might 
provide a method of investigating the nature of millisecond pulsars if these compact stars 
are the central engine of GRBs. 
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Fig. 1. — Evolution of the fireball's Lorentz factor in homogeneous ISM model(a) and wind 
model(b). The solid line is plotted with the model without energy injection, the dotted 
line corresponding to = lO^^G, the dashed line to Bs = lO^^G, the dash-dot line to 
Bs = lO^'^G. The other parameters and values are presented in the text. 
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Fig. 3. — Optical light curve for the R-band afterglow of GRB 000301 C. The data are taken 
from Rhoads & Fruchter^'^'^], Sagar et aL^^^' and Masetti et alJ^^^ and we added 5% systematic 
uncertainty. We fit those data with our homogeneous model: Eq — 1.0 x lO^^ergs, Mq = 5.0 x 
IO-^Mq, Pi = 0.5ms, = 5 X 10W-^ = 0.02, = 5.0 x 10-^ = 3.0 x lO^^G,^ = 3.0 
and Dl — 13Gpc. See the text for more details. 



